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Abstract 

The processes of production and subsequent decay of W- and Z-bosons in e + e" col- 
lisions are considered in a recently proposed modified perturbation theory (PT), based 
on a direct expansion of probabilities instead of amplitudes. In such an approach the 
non-integrable singularities in the phase space, which are intrinsic in the conventional 
PT, appear as singularities in the coupling constant (with subsequent compensation by 
the decay factors of unstable particles) . In the present paper the systematic investiga- 
tion of the modified PT is carried out. The results are compared with the results of the 
conventional approach, based on calculation of the amplitude with Dyson resummation. 
A solution to the problem of the loss of one-loop PT order in the resonance region is 
found. On the basis of this solution the proof of gauge cancellations in any order of the 
modified PT is given. A simple generalization of the fermion-loop scheme is proposed 
which provides a complete description of W-pair production in next-to-leading order 
approximation . 
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1 Introduction 



In many applications of the Standard-Model, connected with the present and future collider 
experiments, one should take into account the effects of the instability of W- and Z-bosons 
(as well as of Higgs boson, top quark etc.) QTJ . In quantum field theory the conventional way 
to take into account an instability consists in the Dyson resummation of the self-energies of 
unstable particles |fjj. This procedure avoids nonintegrable phase-space singularities caused 
by the processes of production and decay of intermediate unstable particles. However, Dyson 
resummation leads to a deviation from the scheme of fixed-order calculations in the frame- 
work of perturbation theory (PT). In gauge theories this results in the violation of the Ward 
identities (WI) and loss of gauge invariance |3|, £|, |[. This fact leads to loss of the control of 
the high-energy behavior of the theory and to the emergence of large errors in the description 
of particular processes. 

In the case of single Z-boson production (LEP1) and within the precision defined by 
one-loop corrections to the vertex functions, the problem of the gauge invariance was solved 
ad hoc (in fact, the mentioned precision implies the approximation of next-to-leading order, 
NLO). The most consistent scheme of calculations was described in [g], where only the gauge- 
invariant contributions to the self-energy were Dyson resummed, while the gauge- dependent 
ones were considered by the conventional PT. As a result, the amplitude could be presented 
as a product of gauge-independent factors, two vertex and one resonant. Nevertheless, this 
result is not universal. Anyway, now it is not clear whether this holds within the next order 
of precision determined by two-loop corrections to the vertex functions. 

In the case of a pair production of unstable particles (LEP2) the amplitude cannot be 
presented in the framework of Dyson resummation in a completely gauge invariant formQ. 
The hopes, nevertheless, for any further progress in such calculations are usually connected 
with the rather general idea of the determination of the minimal set of Feynman diagrams 
that are necessary for compensating the gauge violation by the Dyson resummed self-energies. 
For this purpose within the NLO approximation the fermion-loop scheme was proposed || [| . 
Nevertheless, the bosonic corrections were not taken into consideration in this approach. In 
order to solve this problem a generalization of the fermion-loop scheme was proposed [|10| . 
defined in terms of the formalism of the background-field method. It solved the problem 
of the bosonic corrections and, moreover, it remained in force also beyond the one-loop 
approximation. However, this approach could not solve another problem, which was in 
the fermion-loop scheme, too. This problem concerns the incompleteness of the description 
within the declared precision of description because of the loss of one-loop order of PT in 



the resonance region 11, 12 



Let us consider in more detail the latter phenomenon. The point is that the denominator 
of unstable-particle propagator in the resonance region, p — M 2 = 0(g 2 ), is of order 0(g 2 ), 
but not 0(1). So, in this region the Dyson resummed one- loop self-energy actually makes a 
contribution in the leading-order, but not in NLO. Therefore, in order to complete the NLO 
approximation the two-loop correction to the self-energy should be Dyson-resummed. 

1 It should be noted that there is an alternative approach called the pole scheme, described in || and then 
elaborated in numerous papers. The gauge invariance in this scheme is initially maintained, but, unfortu- 
nately, an algorithm for the evaluation of corrections is not developed completely (see |?| and discussion in 
Sect.|). 
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However, the Dyson resummation of the two-loop correction, made without taking into 
consideration the two-loop corrections to the vertex functions, hardly leaves a chance to 
maintain WL Consideration of all two-loop corrections maybe solves the problem. Never- 
theless, this solution is certainly impractical fLl| , |12 |. Besides, the question remains unclear 



how the two-loop corrections to the vertex functions, which in fact contribute in NNLO, can 
compensate the gauge violation that occurs in NLO. So, anyway the problem of the loss of 
one-loop PT order has to be solved in a practical fashion. 

The present paper proposes the solution to this problem in the framework of the modified 



PT ||13|| . Its basic idea is the expansion of direct probabilities in powers of the coupling 
constant instead of amplitudes. (The amplitudes prior to calculation of the probabilities are 
considered to be full.) Such an order of operations, taken together with ideas of the theory 



of generalized functions [111, allows one to trace the fundamental connection between the 



origin of the phase-space non-integrable singularity and the loss of one-loop PT order in the 
resonance region. (In fact this loss manifests itself as the emergence of an extra singularity in 
the coupling constant instead of the phase-space singularity.) Moreover, while considering the 
probability the contribution of the Dyson resummed two-loop correction (in the amplitude) 
may be reproduced within the given precision in the form of an additive anomalous term. 
Owing to the additivity of this term, it becomes possible to give an independent proof of the 
gauge cancellations in the probability. It is worth mentioning that this result does not mean 
that the inclusion of two-loop corrections only to denominators of "unstable" propagators 
(without the inclusion to vertices) does not lead to violation of WI in the amplitude. This 
means that contributions that violate WI in the amplitude turn out to be beyond the given 
precision in the probability. By its nature this phenomenon is connected with the effect of 
changing of the order of individual contributions in the probability due to the emergence of 
a singularity in the coupling constant instead of the phase-space singularities. 

Below we elaborate the above-stated ideas in any order of the modified PT and give a 
general proof of the gauge cancellations in the probability within the given precision. Notice, 



that the latter outcome was practically anticipated in pioneering work of |13| , presenting 
the modified PT. However the reasoning of this work in the part that concerns the gauge 
cancellations was not complete. Indeed, it was based on a comparison with results of the 
conventional approach, but it overlooked the problem of losing one-loop PT order in the 
resonance region. Besides, it omitted the problem of the difference between remainders of 
the expansions of the amplitude and the corresponding probability. Nevertheless, proceeding 
from only the amplitude it is impossible to estimate in a mathematically correct way the 
remainder of the expansion of the probability, since in the resonance region the expansion of 
amplitude faces the ambiguity of the expression 0/0. 

In order to investigate in detail the above-mentioned problem we first perform the sys- 
tematic study of the modified PT approach. Then we proceed with applying the results in 
the theory of electroweak interactions. In general case we find a practical way to keep a 
fixed precision without violating gauge cancellations in the framework of the background- 
field method. Within NLO we also find the solution in the usual formalism, basing ourselves 
on the results of the fermion-loop scheme. 

This paper is organized as follows. In Sect. |2| a statement of the problem is discussed 
by a simplified example of a single unstable particle production. The basic formulas of the 



modified PT are derived in Sect.^] (as a whole, the content of this section follows [131). The 
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properties of the expansion of the unstable propagator squared are studied in Sect.fJ[ Section 
H discusses briefly the soft-photon problem. In Sect.[| the general proof is given of the gauge 
cancellations in processes mediated by unstable particle production. Section [7| is devoted 
to construction of the generalization of the fermion-loop scheme in NLO approximation. In 
Sect. | the results are discussed. 



2 Unstable propagator in AO: 
statement of the problem 

In this section we consider the structure of the denominator of an unstable propagator. We 
reject, for a moment, all factors in the numerator. Then let us consider the propagator in 
the following form: 

A / \ f 1 

A(a;r) 



M 2 -p 2 -S t - ah(r) - ia/(r) ' w 

Here a = g 2 / (4ir) is the coupling constant squared, r = M 2 —p 2 is the kinematic variable, M 
and p are the mass and momentum of the unstable particle (the mass squared is considered 
without the conventional — iO), £ is the renormalized self-energyQ ah and af are its real 
and imaginary parts with the extracted factor a. Here we do not fix the scheme of UV 
renormalization because the results will not depend on it (see below). By definition, the 
property of instability means that / 7^ in some neighborhood of the point r = 0. Owing 
to causality we assume that / > in this neighborhood. In what follows we assume that 
the size of this neighborhood is of the order of 0(a) and that it contains a solution to the 
equation r — ah(r) = 0. Generally speaking, the function h(r) may be nonzero at r = 0. 

The probability of the process of the production and decay of an unstable particle is 
defined by an integral over some region of the kinematic variable of the propagator squared 
W(a; r), multiplied by some weight function, 

P(a) = J dTip(r) W(a;r), (2) 

W(a; r) ee | A(r) | 2 = 1 (3) 

[t — ah{T)\ + a z j z (r) 

The weight function tp(r) stands for the complementary part of the unitarity diagram with 
respect to the given propagator squared. It includes the photons radiation from the ini- 
tial and final charged states. In the general case ip(r) includes also hardware factors of 
experimental devices. 

By virtue of the property of / 7^ the function W(a; r) is finite and, therefore, is 
integrable in a neighborhood of r = 0. However, in the limit a — > the non-integrable 
singularity 1/r 2 appears in W(a; r). This fact means the impossibility of direct application 



2 In the case of vector bosons the self-energy includes two structures, and p^Pv Formula (1) represents 
the contribution of the first structure only. The effect of presence of the p^p^-structure will be discussed 
latter on. 
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of the conventional PT in the case of processes mediated by unstable particles. Nevertheless, 
the expansion in the coupling constant does exist in the probability P(a). 

Actually, the origin of non-integrable singularity means that the result of integration of 
W(a; t) with a weight function tp(r) includes a singularity in a at a — > 0. If one extracts 
this singularity, then the expansion of the remaining integral becomes possible. In the case 
of a power-like singularity, this expansion will be a Laurent expansion, but not a Taylor one. 
(Let us note that the weight function (f(r) actually depends on the parameter a, including 
it, in particular, as a factor. Therefore, the expansion of the integral ultimately may take a 
Taylor form. However, a priori the kind of the singularity is not known. So in order to study 
the problem, we first consider (p(r) as a function independent from a. We also assume that 
<p(t) is a rather smooth function and generally is non zero at r = 0.) 

Let us study the kind of singularity of the integral. Since the singularity originates from 
integration over small r, we may keep in the functions h(r) and /(r) only their leading terms 
of the asymptotic expansion at r — > 0, i.e. approximate them by ho = h(0) and /q = /(0). 
As a result, up to inessential corrections (which will be calculated later on), we get the 
approximation as follows: 

1 

a h ] 2 + a 2 fl 

By virtue of the homogeneity one can deduce from (4) that the integral of W(a; r) with the 
weight function <p(t) has a singularity of 1/a. Indeed, let us divide the range of integration 
into \t\ > const x a and \r\ < const x a with large enough "const". The integration over the 
first range gives a finite contribution, while the integration over the second range gives the 
above-mentioned singularity (this may be verified by a change of the integration variable). 
Moreover, the coefficient at the singularity is proportional to fo 1 and does not depend on 
ho- In fact, fo may be included into the normalization of a, whereas ho does not make 
contribution in the leading order. Indeed, setting h = does not lead to a singularity in r. 
For similar reasons the weight function <£>(t) gives a contribution to the leading-order term 
SIS db trivial factor y?(0). 

So, despite the fact that the expansion of the integrand is an incorrect operation, the 
expansion of the result of integration is worthwhile. Moreover, some properties of such 
expansion may be determined before the actual calculation of the integral. For systematic 
determination of the properties there is a special method called the asymptotic operation, 
AO |TB|, [IT], . Its key point is the transition to an extended interpretation of the integrand 



W(^)= r. _ U , 2M ,. 2 . (4) 



as a product of a kernel of generalized function on a test function [14, 15]. In fact this means 



that the integral is interpreted as a continuous linear functional on a test function. When the 
integral is well defined (before expanding the integrand) the above-mentioned generalization 
does not lead to any modification. However, after the formal expansion of the integrand the 
new interpretation allows one to give meaning to the non-integrable terms of the expansion. 

Thus, the problem of the expansion of the integrand may be solved basically within the 
method of generalized functions. Next, it is necessary to consider the asymptotic properties 
of the expansion. For this purpose an ambiguity of the extension of the interpretation of 
nonintegrable functions in the integrand is used. Generally, it is well-known (e.g. from ex- 
perience of the UV renormalizations) that elimination of divergences may be accompanied 
by ambiguities. When an integral is determined by the method of generalized functions, the 
ambiguities are described by means of so-called counterterms, which are proportional to the 
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delta-function or its derivatives located strictly in the point of the non-integrable singular- 
ity.0 In the AO framework the coefficients at counterterms are unambiguously fixed by the 
requirement of the reconstruction of the result which should be obtained from an expansion 
of the initial integral. (Let us recall that before expansion of the integrand the integral was 
well defined and there were no ambiguities in it.) Moreover, AO gives a practical recipe 
of the calculation of these coefficients in each order of the expansion prior to a calculation 
of the integral. The resulting counterterms contain complete information about the terms 
which are singular in the parameter of the expansion. Simultaneously the counterterms may 
also contain some non-singular contributions which correct the asymptotic property of the 
expansion. 

In the above example the counterterm that describes the leading term of the asymptotic 
expansion of W(a; r) is cj (afo) X S(t) with c is some numerical factor. In the given case the 
value of c, as well as the appearance of this counterterm, follows from the formula which is 
well-known in the theory of generalized functions, 

a 

lim — — — - = 7T<Mt). (5) 

a^O T 2 + a 2 

Thanks to this formula the expansion of W(ct; r) up to an 0(1) correction is defined un- 
ambiguously and turns out to be the delta-function only. Up to the 0(a) correction the 
expansion is non-trivial. Its most general form is 

W(a; r) = n/(af ) 5{r) + [ 1 /r 2 ] + c S{r) + d(-) 8\r) + 0(a). (6) 

Here [1/t 2 ] is the generalized function defined with some prescription. (The most common 
but not mandatory prescription is the principal value. Later on we omit the square brackets.) 
The last two terms in (6) are counterterms that correct the contributions of the first two 
terms. 

In the general case the complete determination of the generalized function 1/t 2 may be 
done by the following way ||15|| . First, one makes two subtractions in the test function (p(r) 
in some neighborhood of r = by replacing <p(r) to (f(r) — y?(0) — rip'(0). As a result, 
the non-integrable singularity in 1/r 2 becomes compensated. (In fact this is one of possible 
prescriptions.) Then in order to describe the ambiguity emerging with this subtraction, one 
has to add two counterterms to 1/t 2 . One counterterm is proportional to the delta-function 
and the other one is proportional to its first derivative (both counterterms correspond to 
the above subtractions). The coefficients at the counterterms must be determined in such a 
way |1| [L7|, [18|| as to guarantee the asymptotic properties of the expansion of the integral at 
order 0(1). Their values depend on h and /, and on the choice of the prescription in 1/r 2 , 
but the sum of all terms will not depend on the prescription. 

The above process may be continued. The next term of the formal expansion of W(a; r) 
is 2ah(r) x 1/r 3 . For its complete determination one needs three counterterms which include 
the delta-function, its first and its second derivatives. The coefficients at counterterms are 
fixed by the requirement of the asymptotic properties of the expansion. The practical recipe 
of the calculation of the coefficients is presented in the next section. 

3 Let us emphasize that the introduction of counterterms is a general place in the theory of generalized 
functions (see e.g. [|l5)). Actually this idea was used by N.N.Bogoliubov p(| for establishing the R- 
operation. In the AO context the term "counterterms" was introduced |l6|. in order to emphasize the 
analogy with the theory of UV renormalizations. 
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3 Calculation of counterterms 



Let us show the technique of the calculation of counterterms on the example of the AO 
expansion of W(a;r) up to 0(a 2 ). Since the leading term of this expansion is of order 
0(a~ l ), the mentioned precision is sufficient for the determination of the next-to- next-to- 
leading order (NNLO) approximation. Such a precision seems to be sufficient for most of 
the practical applications. It is sufficient also for understanding the properties of the AO 
expansion of W(a; r). 

The general structure of the AO expansion of W(a; r) within the considered precision is 
as follows: 

W(a;r) = \ + ^p-a + E(r) + 0(a 2 ) . (7) 

Here the first two terms represent the result of the formal expansion of W(a; r). For defi- 
niteness, we treat the poles with respect to r in the sense of principal value. Below we recall 
two equivalent definitions of the principal value: 



V7>i- = i 

T n 2 



(r+iO) n (r-iO) 



ya-l d n 

In r 



{n-l)\ dr n 



The derivatives in the last expression are understood in the sense of generalized functions, 
i.e. they must be switched to the test function via formal integration by parts. 

The quantity E{r) in (7) represents the sum of counterterms. In this case they are 
proportional to the delta-function, its first and its second derivatives: 

E{r)= £L±*i *(»)(,-). (9) 

n = U - 

In what follows we assume that h(r) and f{r) together with their second derivatives are 
regular functions in some neighborhood of r = O.f] At first, for simplicity we suppose that 
the functions h and / include one-loop contributions only. A generalization to the multi-loop 
case is considered below in this section. 

We wish to define a procedure for the determination of the coefficients c n , n = 0,1,2, 
for an arbitrary (p(r) function, which decreases rapidly enough at infinity. This procedure 
should provide the following: 



J dry?(r) 



T z T 6 



0(a 2 ). (10) 



Note, that for the solution of this problem we do not have to know all information about 
the function h(r) in the third term in square brackets. We need only have knowledge about 
three terms of its asymptotic expansion, 

h(r) = h + rti + (r 2 /2) h% + o(r 2 ). (11) 



4 If the unstable particle is able to interact with massless particles (photons) , this property does not occur. 
Nevertheless, supplying the photons with mass we restore the analyticity inside a neighborhood defined by 
the generated mass gap. This is enough for our purposes (see discussion in Sect. ||). 
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This property follows from the fact that the remainder o(r 2 ) cancels the nonintegrable sin- 
gularity 1/r 3 in (10). 

Now let us substitute (9) into (10) and, following |16[| , present the test function ip(r) as a 
linear combination of three basis functions Lp n {r), n = 0,1, 2, such that <pVf\0) = 5%, where 
5* is the Kronecker symbol. As a result we get 



-f-UHJ 

z n = J dry? n (r) 



W(a;r)- 



2h(r) 



a: 



+ 0(a 2 ). 



(12) 



It is obvious from (12) that within the given precision the coefficients c n do not depend on the 
choice of test functions. Indeed, any other test functions (p n (r) possessing the same property, 
V?n (0) = ^ni ^ ea, d to the coefficients c n instead of c n . However, the difference c n — c n amounts 
to a quantity of order 0(a 2 ), because the difference between the corresponding integrals is 
controlled by the test function (p n {j) —ip n {r), which equals zero at r = together with its first 
and second derivatives. It is easy to see that the integral in formula (12) with such a weight 
function is of order 0(a 2 ). Due to this property, without loss of generality one may choose 
the test functions to be the step-like ones. Namely, we may set ^(t) = r" x &(\t\ < A). 



Then, with n 



0,1,2, 



we have 



+A 



drr 11 



-A 



W(a;r) 



1 2h(r) 



-a 



+ 0(a 2 ). 



(13) 



Formula (13) basically solves the problem stated above. However, it is still too compli- 
cated for practical usage since, generally speaking, through W(a; r) it contains the depen- 
dence on the unknowns functions h{r) and fij). Moreover, (13) contains a lot of superfluous 
information because the integral in r.h.s. includes contributions beyond the given precision. 
In particular, the dependence on the cutoff parameter A is of this kind. 

i.e. by the specific trans- 



The problem may be solved by homogenization procedure [IS 



formations in the integrand. Namely, let us do the substitutions r — > £r, a — > £a. Then the 
result we expand in powers of £, and in the end we set £ = 1. Each term of this (secondary) 
expansion is proved to be a homogeneous function of r and a. So it leads to a strictly def- 
inite contribution in the power in a to the integral, which now may be considered without 
the cutoff. The first term of this expansion gives the leading-order contribution to the 
coefficient c n , 

1 1 2/ k) 



drr n 



(r - ah ) 2 + u 2 fl 



-a 



(14) 



The next term of the expansion of homogenization gives the correction term c£\ etc. On 
counting of powers we have ~ a n_1 , ~ a n , etc. Adding the necessary number of c$ 
we get c n with the required precision. 

Let us emphasize, that the integral in (14) is convergent at infinity. The convergence 
takes place also for other c$. Moreover, this property is a general result of the application 
of the homogenization [I8|. It is worth recalling that the singular terms 1/r 2 and 1/r 3 
in (14) are defined in the sense of principal value. Nevertheless, the prescription may be 
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changed (in all above- derived formulas). As a result the coefficients c$ will change, but the 
asymptotic properties of the expansion (7) will be conserved. 

After the corresponding calculations we come to the following result (for the first time 
obtained in p3|): 



Co 



7T ir(h' f - hofo. 



a/o fo 



+ 73 

Jo 



nh n (2h h' f - ftg/' Q - fjf ) tt (§ - / 2 ) 

ci = —r- H 72 a ' c 2 = 7 a • I 15 ) 

Jo /o Jo 



Here the subscript means that the corresponding quantities are defined at r = 0, while the 
superscript means we are taking derivatives. For example, h' Q = dh(r) / dr\ T=0 , etc. 

The derived result may be written in a more compact form if the quantities c n in (9) are 
understood as functions of r. In this case, instead of (15), one can obtain 

7r(h 2 -f 2 ) 

c = — r , ci = — , c 2 = a . (16) 

af f f 

Here Co i,2, h and / are functions of the r. The equivalence of these two forms, (15) and 
(16), follows from the relations /(r) 5'(r) = f S'(r) - /^(r) and f(r)8"(r) = f 8"{r) - 

V'o&'{T) + f»&{T) . 

The above-derived results may easily be extended to the case of multi-loop contributions 
in h(r) and f(r). In this case one should perform the conventional expansion in a in formulas 
(15) or (16), in which h and / are understood with the presence of multi-loop contributions 
p~3|] . The simplest way to prove this statement is as follows: assuming that h and / are 
the full functions we repeat all the same reasoning as was done above, except that during 
homogenization we modify the scaling by setting a n — > £a™ in the n-loop contributions 
(instead of a n — > £ n a n ). As a result the higher-loop contributions become identified with 
the one-loop ones, and formulas (15) and (16) are restored automatically. 



4 The properties of AO expansion of W(a; r) 

Let us discuss now the properties of the AO expansion of W(a; r). At first we are interested 
in the characteristic which will be useful in proving gauge cancellations in the probability. 
It should be emphasized that a solution to this problem is not obvious a priori due to the 
specificity of the Feynman rules in the modified PT approach (see formulas (7), (9) and 
(15)). Then we will examine the self-consistency properties of the expansion, such as a non- 



sensitivity of the results to the ambiguities in the definition of the gauge-bosons masses JsT 
and consistency with the requirements of UV renormalization. 

We start from a rather methodological problem consisting in the explicit demonstration 
of the independence of the results of an AO expansion from the sequence of the expansion 
in the higher-loops. In other words, we show that the result of the AO expansion of W(a;r) 
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will be the same if instead of (1) one will start from the following formula of the incomplete 
Dyson resummation: 

A( t) = 1 = 1 a 2 E 2 + a 3 £ 3 

Here a n S n (r) is the n-loop contribution to self-energy. By squaring (17) we get the incom- 
plete expansion for W(a; r), each term of which is integrable in the conventional sense: 

W(a; t) = Wi(a; r) + [(a 2 £ 2 + a 3 S 3 ) W u (a; r) + a 4 (£ 2 ) 2 W 12 (a; r) + h.c." 

(18) 

+ a 4 |E 2 | 2 W 2 (a;r) + 0(a 2 ). 

Here we have introduced three new functions: W n (a; r) = W^a; r)Ai(a; r), W 12 (a; r) = 
W^a; r)A 2 (a; r), and W 2 (a;;r) = W ^a; t)W ^a; r) , where W^a; r) and A 1 (a;r) are the 
same functions as W(a; r) and A(a; r), but with only one-loop self-energies Dyson resummed. 
By considering these new functions as generalized functions, each term in (18) may be in turn 
completely AO expanded. Repeating the reasonings of Sect.|], one can show that the leading 
term of AO expansion of W n (a; r) has the behavior of 1/a 2 . So, after integration with the 
weight function (f(r) the first term in the square brackets in (18) makes a contribution of order 
0(1). The leading terms in Wi 2 and W 2 both are 0(l/o; 3 ). So, the second term in the square 
brackets and the last term in (18) are of order 0(a). By a similar reasoning one can show 
that the neglected terms in (18) are of the order 0(a 2 ). So in view of W^a; r) = 0(a _1 ), 
formula (18) approximates W(a; r) within the NNLO precision. 

It should be noted, that the above reasoning is valid only as long as the functions £;(t), 
/ = 2, 3, ■ • •, and several of their derivatives are regular functions in some neighborhood of 
t = 0. If this is not the case, then the products of X/(t) on must be considered as 
new generalized functions, of which the properties should additionally be investigated. Such 
situation occurs when the unstable particle interacts with massless particles (photons). This 
difficulty may be eliminated by the inclusion of the regularizing mass for massless particles, 
because then functions S;(r) become regular in a vicinity of r = 0. 

Assuming this trick let us now consider the complete AO expansion of W^. For brevity 
we omit the corresponding derivation since it is similar to the one considered in the previous 
section. The desirable accuracy of the expansion is controlled by formula (18), so that Wn 
should be expanded up to 0(1) corrections, and W i2 and W 2 up to 0(aT 2 ) corrections. 
However, in further calculations major precision is required. So, let us consider 

W u (a;r) = E(r) + \ + 0(a), W 12 (a;r) = E{r) +0(1), W 2 (a; r) = E(r) + 0(1). (19) 

T 

Here in all three cases E(r) is defined by (9), but with different coefficients c n . In the 
compact form, when c n are defined as functions on r, we get: 

Wn - Co "2^p' Cl ~ laf ' C2_ W ' (20) 

W ■ r ~ * r - r - + f - h 2 ) _ 
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Notice here that the coefficients singular in a do not depend on the prescription for the 
poles in r, since in the examples considered above the non-integrable terms (for which the 
prescription is needed) are non-singular in a. 

Based on the derived results we formulate the following properties of the expansion. (All 
of them can be verified by direct complete expanding and comparing the results.) 

Property 1. The incomplete expansion (18) is equivalent within the given precision to the 
complete AO expansion of W(a; r). 

The next properties are non-trivial and some additional argument should be given. 

Property 2. Within the given precision an incomplete expansion of W(a; r) is equivalent to 
its complete AO expansion if this incomplete expansion includes as Dyson resummed all con- 
tribution non-zero at r = to ImE^r). All other contributions to Ei(t) may be transferred 
from denominators to numerators in the sense of a conventional expansion. Moreover, only 
a finite number of terms of the latter expansion are relevant within the given precision (see 
remark at the end of the proof). 

We perform the proof in two steps. At first we will show that all contributions to Ei(r) 
zero at r = without loss of precision may be expanded in the conventional sense. Then we 
will show that the same operation may be done also for the whole of the real part of Ei(r). 

So, let £i(t) = Soi(r) + Ei(r), where by definition Ei(0) = 0, but E O i(0) + 0. Since 
in some neighborhood of r = Si(r) is a correction to E 01 (r), its contribution may be 
expanded like E n (r) with n > 1 in (18): 



W(a;r) =Wi(a;r) + 



(«S! + a 2 S 2 + a 3 S 3 ) Wn(a;r) 
+ (a 2 E 2 + 2a 3 SxS 2 + a 4 E 2 ) W 12 (a; r) + h.c. 



(23) 



+ 



a 2 



Si " + a 3 Ei E 2 +h.c. + a 4 |E 



J 2 



W 2 (a;r) + 0(a 2 ) . 



Here symbol "*" means complex conjugation. In the denominators of only the E i(t) * s 
Dyson resummed. The remainder in (23) is estimated in the AO sense. The AO expansions 
of Wi, Wn, W12 and W 2 within the required precision were previously described in this 
section. 

The proof of (23) may be given noticing that within the given precision the quantity 
ckEi gives a non-zero contribution only if it is raised to not more than the second power. 
In fact, against the background of the regular terms this property is obvious. If E 2 (or 
|Ei| 2 ) is multiplied by a counterterm originating from W^, the result is non-zero only in 
the presence of second- or higher-order derivatives of the delta-function in the counterterm 
(otherwise there acts the property E^O) = 0). Consider the functions = [Wi] m A™, 
n > 0, m > 0, n + m — 1 is the number of self-energy insertions on one side of the cut 
of the diagram of unitarity, m — 1 is the number on the other side. In these functions 
such counterterms appear in order a^( n + 2m ^ 1 ) x a 2 , and only those functions are to 
be multiplied by E 2 (or |Ei| 2 ) which satisfy the condition n + 2m — 2 > 2. Note that in 
both insertions of E 1; each gives a factor a. Of the other remaining n + 2m — 4 insertions 
of Efe, k > 2, each gives a factor not less than a 2 . Thus all considered terms result in 0(a) 
contribution. While extending the above reasoning to the third and the higher powers of 
aEi, we easily see that they give nonzero contributions of 0(a 2 ) only. 
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The above result may be generalized to the real part of £i(t). In this case £i(t) in (23) 
is to be defined by ImS^O) = with ImS i(0) 7^ 0. The basis for this generalization is 
the observation that the real part of the self-energy does not contribute to the leading-order 
term of the AO expansion of W(a; r). Let us remark, that at first sight formula (23) with 
this modification should look much more complicated, because from a formal point of view 
it should contain an infinite series of terms of [a ReSi(r)]™ x Wi„, each of which is of order 
0(a _1 ), etc. However, forming the groups with other functions all superfluous terms 
must mutually cancel. The mentioned groups will be formed by virtue of relations of the 
type 2ReW 12 + W 2 1 = 0{a~ l ) [not 0(a~% etc. 

It is worth noticing that the above result permits one to expand the bosonic corrections 
to self-energy of W- and Z-bosons within the finite number of terms. (This is because 
the mentioned bosonic corrections possess the property lHi£k os (o) = 00 P3-) This is a 
nontrivial result since according to the formal power counting all these terms are of the same 
order in a and, consequently, should be explicitly taken into consideration. 

Property 3. There is the following approximation of W(a; r) up to 0(a n ) corrections: 

W(a; r) = W [n] (a; r) - a 1 - 1 ^] n6(r) + 0(a n ). (24) 

[ImSi(0)J 

Here W[ n ](a;r) stands for the propagator squared with Dyson resummed contributions to 
the self-energy up to n-loops. The second term includes the (n + l)-loop correction which 
is necessary (with the appropriate factors) within the indicated precision. Formula (24) 
follows immediately from the obvious generalization of (18) to the case of any n taking into 
consideration the first result in (20). 

It is worth noticing that the second term in (24) makes a contribution in the highest order 
within the given precision. Consequently, any 0(a) -variation of the mass is worthless in this 
term, because the effect may be referred to the neglected term of order 0(a n ). Actually 
this is a very important observation, since in the case of electroweak theory it permits one 



to leave out of account the ambiguities in the definition of the gauge-bosons masses |2T|. In 
particular, one may consider the factor ImE n+1 (0) to be defined at the complex- value pole 
position in the spirit of the pole-scheme (as a result ImS„ + i(0) becomes completely gauge- 
invariant), or simply as the n-loop correction to the width of the unstable particle divided 
by its mass. 

Formula (24) represents the quantitative characteristic of the property of the loss of one- 
loop PT order in the resonance region. Indeed, as long as Wu(a; r) ~ a -1 at a — > 0, the 
second term in (24), represents the n-th order correction, but not the (ri + l)-th one, as might 
naively be expected. Since the second term in (24) cannot be obtained from the analysis of 
an amplitude, it is pertinent to call it the anomalous additive term. 

Property 4- The above expansions possess the following transformation property when the 
argument r is shifted by a quantity of order 0(a): 

W(a; t) = W(a; r - am 2 ) h(r _ am 2 ) _ ft(T) _ m 2 . (25) 

/(r- Q m 2 H/(T) 

Here W(a; r) stands for any incomplete or complete AO expansion of the considered above 
functions; the quantity m 2 is of order 0(1). 
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Property (25) means the non-sensitivity of the entire formalism with respect to variation 
of the mass shell within 0(a). Moreover, it means independence of the formalism from the 
UV renormalization scheme. Indeed, at the one-loop level the transition, for example, from 
the MS scheme to the on-mass-shell (OMS) scheme is described by f] 

S OM5 (p 2 ) = £(p 2 ) - Re£(M 2 MS ) - (p 2 - M 2 OMS ) x Re£'(M 2 MS ) , 

M 2 oms = M 2 - ReS(M 2 MS ) . 26 

It is obvious that (26) belongs to the class of transformations covered by (25). The trans- 
formation at the multi-loop level is controlled by formula (18). So this can be done in 
accordance with the standard recipes of the UV renormalization, which do not depend on 



presence of the "infrared" counterterms located on the mass shell (see also [17] and references 
therein). Let us note that the transformation to another scheme of UV renormalization may 
be done (speculatively) before the squaring of the amplitude and the AO expansion. The 
consequent AO expansion by no means "feels" in what scheme the Green functions have 
been renormalized. 

The property (25) is trivial in cases of non-expanded W or its formal expansions. The 
nontrivial aspect is that it remains valid for the counterterms E(t), too. However, this also 
can be understood if one notes that the transformation r — > r — am 2 does not affect the 
structure of the homogenization at the scaling r — > £r, a — > £a (see Sect.0). As a result, 
the property (25) is valid in the most general case. 

Finally, we present one more property which represents doubtless independent interest. 
Property 5. From (25) follows the next recurrent formula for the coefficients c n , 



1 

Cn-l — — 

n 



1 d ^ n /,( r +i) d {r+ i) d \ d 



adho r e V dht ] u dfi r) ) dM 2 



(27) 



Here c n = c n (M 2 ; a; ho, ■ ■ ■ tiQ N ~ n ' 1 ; /o, ■ ■ ■ fo N ~ n ^) are the constants independent from r, n 
runs values < n < N, where N is the maximum degree of the derivative of the delta- 
functions in E(t). Formula (27) is written down with possible dependence of the coefficients 
c n on the parameter M 2 . An essential point for the derivation of (27) is the expansion in 
powers of a of the delta-function 8{r — am 2 ) and its derivatives in r.h.s. of (25). 

The practical value of (27) consists in the opportunity to check the results of calcula- 
tion of counterterms, or to determine the "lower" coefficient c n _i to within 0(a L ) if the 
"higher" coefficient c n is known to within 0(a L+1 ). 



5 Massless particles exchange 

The problem of taking into consideration massless particles (photons) requires a special 
analysis because massless-particle contributions to the self-energy of a massive (unstable) 
particle involve the singularity of rx ln(r — iO). The first derivative of this expression is not 
defined at zero. So already the first corrections to coefficients (15) become uncertain. 

5 Remember, we do not consider contributions to the numerators of propagators. Therefore, we disregard 
the multiplicative wave function renormalization. 
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One way to solve this problem is to introduce a regularizing mass for massless particles (a 
soft-photon mass). Then the non-analyticity at r = disappears, and this fact opens a way 
to use without problems the above-derived formulas of the AO expansions. In a properly 
defined probability, taking into account radiation of the real soft photons, all contributions 
singular in the photon mass should cancel among themselves. That leads to continuity of 
probability with respect to the photon mass. So, while solving qualitative problems, one 
could not worry about the presence of the photon mass, because in the final results the 
dependence on it may be eliminated by the ordinary passage to the zero- mass limit. 

It should be stressed that in the above discussion the photon mass should not necessarily 
be infinitesimal. Moreover, it may be chosen to be so large as far as is needed for the solution 
of a particular qualitative problem. In what follows we use the photon-mass insertion for the 
consistent determination of the orders of expansion in the coupling constant of the photon 
contributions to the probability. After completion of all cancellations and after the passage 
to the zero-mass limit these orders will not change. It is worth noticing that the photon 
mass may be introduced in such a way as to guarantee validity of WI. This can be seen by 
consideration of the problem in the Stueckelberg formalism]^ Moreover, a photon mass (and 
a gluon mass, as well) may be introduced in a totally gauge-invariant fashion ||23|| . So, the 
photon mass will not brake down gauge cancellations. 

Another way |13|] to solve the problem of massless particles is based on usage of the 
regularization property of the parameter a. This method is deeply involved in the context 
of the AO expansion. It should be noted, however, that it is able to cure only those IR 
divergences of which the origin is connected with the emergence of an additional singularity in 
a at a — > 0. In reality such divergences appear only in the diagrams where the soft momenta 
of massless particles come into the propagator of unstable particle considered on shell. Let 
us note that cancellation of these IR divergences means cancellation of the corresponding 
singularities in the coupling constant, and vice versa. The idea of the method of consists 
of a stepwise expansion of the full Green functions squared: first in the contributions of 
the massless particles only, and then in other vertices using the AO technique. (The first 
step will not lead to an infinite series of equal orders in a in the resonance region owing to 
Property 2 and the property ImS5 >os (0) = for the massless-particle contributions.) 



Acting in the framework of the method of [fL3|l , after the first step one gets the modi- 
fied Green functions with unstable propagators not containing contributions of the massless 
particles. However, this simplification is not indisputable, since there will appear a lot of 
configurations with the explicit contributions of photons for which some special counterterms 
will be needed. (In fact they will regularize the products S; x W^; see the discussion in 
the previous section.) Basing ourselves on considerations of unitarity one can justify can- 
cellation of the considering class of IR divergences [[13] . However, the property of the gauge 
cancellations still remain to be proved. So this method of handling the IR divergences seems 
not too good, especially since it does not provide a complete solution of the problem and 
requires a lot of additional efforts. 



The author is grateful to A.A.Slavnov for indication of this fact 
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6 The gauge cancellations 



Now we are ready to give a proof of the gauge cancellations in any order of the modified 
PT. We present our argumentation in a rather general form, applied basically to any un- 
stable particle in electroweak theory. The main idea is to separate in the probability the 
contributions certainly possessing the property of gauge cancellations from the problematic 
contributions. That allows us to concentrate on study the problematic contributions only. 

Let us begin with some preliminary notes. First of all we determine the photon-mass 
regularization for IR divergences. The advantages of this method have been discussed in 
the previous section. Here we stress that this method permits one to consider soft photons 
like ordinary massive particles. That means that their contributions either may be referred 
completely to vertex blocks, or, in other cases, they suppress the resonant behaviour of the 
corresponding contributions. Really, an exchange by a massive particle suppresses the on- 
shell propagating of unstable particles both before and after the emission of the massive 
particle. As a result, the problem of instability is reduced solely to configurations which 
include only the pairs of unstable propagators of equal mass and momenta, placed on both 
sides of the cut of unitarity. It is worth remembering that the mass insertion for photons 
will not violate gauge cancellations (see previous section). Moreover, since the probability 
is a continuous function of the photon mass the dependence on it may be eliminated in the 
final results by the ordinary passage to the zero- mass limit. 

The next note concerns the unphysical pole contribution to the vector boson propagators. 
Assuming the parameterization for the self-energy 

— Hg tlv + S L p^p u , (28) 

we come to the following explicit expression for the full propagator in i?^-gauge: 

Here the first term represents the product of the spin factor on the scalar propagator A (a; r) 
introduced in (1). The second term describes the unphysical contribution. In the conven- 
tional PT it is cancelled due to WI by contributions of the other unphysical states. Con- 
sequently, the second term in (29) does not lead to non-integrable singularity in the phase 
space. In the modified PT this fact allows us to take into account the second term in the 
conventionally expanded form. Then, the property of gauge cancellations will be controlled 
by the properties of the first term, taken in an absolute value and having been squared. 
Note that the cross terms resulting from the squaring of the propagator (29) will not lead 
to non-integrable singularities in the phase space. So they also may be taken into account 
as conventionally expanded. 

Now let us show that the first term in (29), being taken into consideration in framework 
of the modified PT, does not break the gauge cancellations in the probability. Remember 
that in accordance with (24) its contribution may be presented in the form of a sum of two 
terms, where one is conventional and the other is anomalous. Let us start from analysis 
of the anomalous term contribution. Due to its additivity the simplest way to perform the 
analysis is to assume for a moment that the squared unstable propagator consists in this 
anomalous term only. 
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Firstly we consider the case with a single unstable particle production. Let us note that 
the complementary part of diagram of unitarity with respect to the given propagator squared, 
by virtue of the delta-function without a derivative in the anomalous term, is taken strictly 



on mass shell. Moreover, the details [21 1 of the definition of the mass shell are inessential 



here and may be omitted. (This is because these details have meaning of corrections, but 
the anomalous term describes the highest order within the given precision.) From these two 
facts follows the conclusion that the sum of the complementary parts is gauge invariant, since 
actually it coincides (up to factors) with the product of two S'-matrix elements squared, with 
one describing the on-shell production of the unstable particle and another one describing 
its decay. Furthermore, in the anomalous term itself one may neglect the gauge- dependent 
contributions, if in this case they are in ImS n+1 (0), since they are beyond the given precision 
(see Property 3 and the discussion therein). So the total sum of unitarity diagrams that 
include the anomalous term within the given precision is gauge invariant. 

In the case of multiple production of unstable particles let us single out any unstable 
particle and, again, consider its contribution of an anomalous term. Since the anomalous 
term is of the highest order, the contributions of all other unstable particles should be 
considered in the leading order of the AO expansion. By virtue of (7), (9), and (15) the 
leading-order contributions are determined by the delta-function without a derivative with 
a coefficient including the one-loop fo only as a non-trivial factor. In view of the gauge 
invariance of the one-loop fo (see, for instance, ||) follows, once again, the gauge invariance 
of the total sum of complementary parts of the diagrams of unitarity. So, the gauge invariance 
in this case occurs, too. Due to additivity, a similar reasoning may be repeated for any other 
unstable particle of multiple production, and in each case we will reach the property of gauge 
invariance of the corresponding contribution to the probability. 

Now let us proceed to the analysis of all other contributions, which do not at all include 
an anomalous term. Such contributions are controlled by the first term in (24). Also, all non- 
resonant contributions belong to this class of contributions. As a whole, these contributions 
are described by the amplitude squared determined in the conventional approach with Dyson 
resummation. The property of gauge cancellations in the amplitude so determined has al- 
ready been shown |10[ in the framework of the background-field formalism. Let us emphasize 
that the only condition which has to be observed when considering these contributions is the 
Dyson resummation of all corrections to the self-energies up to n-loops together with taking 
into consideration all n-loop corrections to the vertices. So, treating W[ n ](a;r) in (24) in 
such a manner, we automatically get the property of gauge cancellations in the probability^ 

The above discussion completes a construction which provides both the gauge cancel- 
lations and the necessary precision of the description in the sense of an expansion in the 
coupling constant. Basing ourselves on this result one may do the next step: proceeding 
to the complete AO expansion of the probability. Then the property of gauge cancellations 
within the given precision must take place as well, because AO only expands in powers of 
the coupling constant the expression which possesses the property of gauge cancellations. 



It should be noted that within the background-field formalism after all gauge cancellations there may 
remain some residual dependence on the quantum gauge parameter. This phenomenon has been considered 
in the original work of friof , where it has been stressed also that this residual dependence will not affect 
the high-energy behaviour of the amplitude. An idea for how to solve the problem of the residual gauge 
dependence is discussed in SectR 



16 



With this note we complete the general proof of gauge cancellations in the approach of the 
modified PT. 



7 Generalization of the fermion-loop scheme 

If the analysis may be restricted by NLO precision, for example for the case of W-pair pro- 
duction studied at LEP2, the gauge cancellations may be proved within the usual formalism, 
without applying the background-field method. The key point is the well-known result on 
gauge cancellations in the so-called fermion-loop scheme 0, Recall that it consists of 
including all fermionic one-loop corrections in tree-level amplitudes and Dyson resumming 
self-energies. The difficulty of this scheme is twofold. First, it is not known how to incorpo- 
rate the one-loop bosonic corrections into this scheme without spoiling the gauge invariance. 
Second, there is a problem with gauge invariance while taking into account the two-loop 
corrections to self-energies in the denominators of unstable propagators, which is necessary 
in the resonance region. 

Both these problems may be solved within the modified PT approach with usage of the 
AO technique. Let us begin with the problem of the two-loop corrections. As we have seen 
above, they can be taken into consideration without breaking gauge invariance by adding 
the anomalous terms to the probability. In view of (24) that may be done by means of the 
formula 

W(a; r) = W {1] (a; r) - ^(^ 2 ^ + ° (a) ' (3 ° } 



Here W[i](a;r) represents the unstable propagator squared with Dyson resummed all one- 
loop corrections. Remember that Wm(a;; r) = 0(a~ 1 ) at a — > 0. Therefore, the second term 
in (30) describes the highest-order (NLO) correction. So, any ambiguities in its definition, 
in particular those which concern the problem of the gauge invariance, may be referred to 
the discarded terms of order 0(a). (See Property 3 and the discussion therein.) 

The problem of one-loop bosonic corrections may be solved, too. Let us group these 
corrections into two classes. To the first class we refer the corrections to self-energy of 
unstable particles. To the second class we refer the corrections to the vertex factors, the 
corrections due to exchange processes, and due to the real (soft) photons. 

The corrections of the first class can easily be taken into consideration by means of the 
fact that the imaginary parts of the on-shell bosonic corrections to the self-energies of W- 
and Z-bosons are zero ]5|, 22 1. Owing to this fact and Property 2 of Sect. ^ these correc- 



tions may be transferred from the denominators of unstable propagators to their numerators 
without loss of precision. Moreover, in an OMS-like scheme of UV renormalization, where 
the renormalized self-energies satisfy the condition ReEi(O) = ReS^O) = 0, one has the 
following approximation: 

W [ i](a;r)='WiF(a;T) + 0(a). (31) 

Here W 1 ^(«;r) represents the propagator squared where only the fermionic one-loop cor- 
rections are Dyson resummed. Substituting (31) into (30) we come to the formula which 
leads to gauge cancellations. In order to see this, we should only repeat the reasoning of the 
previous section keeping in mind the gauge invariance in the fermion-loop scheme. 
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However, the bosonic corrections of the second class have not yet been taken into account. 
In order to do that let us make use the fact that within the NLO approximation and in the 
presence of the corrections, the quantity Wlf(oi;t) must be taken into consideration in 
the leading-order approximation only. Remember that in this approximation Wip(a;r) = 
7t/(q;/of) x S(r), and this expression is explicitly gauge invariant. Moreover, since the delta- 
function enters without derivatives, the sum of all factors appearing in its presence in the 
unitarity diagrams are also gauge invariant (see previous section). 

Thus, we come to the following recipe of the generalization. We formulate it having in 
mind the total cross-section for the typical LEP2 processes CC10, CC11 and CC20, which 
have been studied in the framework of the fermion-loop scheme 0. In fact, the generalization 
consists of adding to the probability two types of corrections. 

The corrections of the first type describe the anomalous contributions. In the NLO 
approximation they look like the cross-section of the pair on-shell production of unstable 
particles taken in the leading-order approximation, times the leading-order decay blocks of 
unstable particles, and times the "anomalous" factor. Indeed, the presence of the anomalous 
factor means that all other factors should be taken in the leading order. However, in the 
leading-order approximation only the double- resonant subprocesses contribute to e + e~ — > 4/. 
(In fact these subprocesses are of CC03 class.) That follows from the fact that only these 
subprocesses include the factor 1/a 2 originating from the product of two unstable propa- 
gators squared. Now, let us remember that the additive anomalous term in (30) includes 
the delta-function as well as Wi^(a;r) does in the leading-order approximation. So, the 
anomalous terms always contribute "on-shell" . 

The corrections of the second type may be represented, again, as the cross-section of the 
pair on-shell production of unstable particles, multiplied by their decay blocks. However, 
instead of "anomalous" factors, they include the bosonic corrections to the vertex blocks 
and the corrections due to the real soft-photons. The examples of the unitarity diagrams 
which generate contributions of this type are shown in Fig. |l[ With the non-zero mass of the 
photons all these diagrams include exactly two pairs of unstable propagators of identical mass 
and momenta in both sides of the cut. (The only exception, configuration (l.e), is discussed 
below.) Therefore, they include two delta-functions of the corresponding kinematic variables, 
which make these configurations on shell. The sum of all such configurations represents 
the product of the cross-section of the pair on-shell production of unstable particles and 
their decay blocks. Since these quantities are continuous functions of the photon-mass, the 
dependence on it may be eliminated in the very end of the calculation by the ordinary 
limiting procedure. Notice that the mentioned property of continuity follows directly from 
the well-known theorem for cross-sections with the real soft-photon contributions. 

Now let us discuss the mentioned above exceptional configuration (l.e). Strictly speak- 
ing it should not be considered to belong to the unitarity diagrams, since it describes the 
self-energy correction to the unstable propagator, which has already been taken into consid- 
eration by formula (31). Nevertheless, while considering the cross-section of the pair on-shell 
production, one has to take into account the virtual soft-photon insertion to the external 
legs, which is due to the wave function renormalization. The configuration (l.e) was added 
to the list of diagrams of Fig. |T] only in order to indicate this fact. 

It should be noted that among the diagrams of Fig. [I] there are both factorizable and non- 
factorizable configurations in the sense of the classification of |7], O, |T2| . Nevertheless, at the 
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Figure 1: Examples for factorizable and (quasi) non-factorizable corrections to W-pair pro- 
duction in e + e~ — > 4/(7) which have to be taken into consideration in NLO approximation. 
(The dashed lines denote massive bosons or soft-photons. The continuous thin and thick 
lines represent the initial/final fermions and the unstable W-bosons, respectively. The verti- 
cal dot-dashed lines indicate the cut of the diagrams of unitarity. The hatched areas denote 
the lowest-order Green functions for the production of the virtual W^-boson pair.) 



(a) (c) (e) 




(b) (d) 



Figure 2: Examples for non-factorizable corrections to W-pair production which contribute 
beyond the NLO approximation. 
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intermediate step, when the soft-photons are supplied with the mass, the non-factorizable 
configurations of Fig. [I] become factorizable. All other non-factorizable corrections, which 
do not provide this property, produce such configurations that include less than two pairs 
of unstable propagators of identical mass and momenta in both sides of the cut. Therefore, 
they do not include the leading-order factor 1/a 2 . In view of the presence of the additional 
factor a due to the bosonic corrections they contribute beyond the NLO approximation. The 
examples of configurations of this type are shown in Fig. |2|. It is worth noticing that they 
may not be considered as corrections to the cross-section discussed above of the pair on-shell 
production of unstable particles, or to their decay blocks. 

The above discussion leads us to the following formulas for the total cross-section: 

i 

a(s) = / dz <j>{z\ s) <j(zs), (32) 



o 

8 (v^-v 7 ^) 2 

<t(s) = J ds + J ds_ <7 (s ; s+, s_). (33) 
o o 

Here a(s) is the experimentally measured cross-section at the center-of-mass energy squared 
s. (p(z; s) is the "flux" function?] describing the contributions of the initial- and final-state 
photon radiations with large IR and collinear logarithms, a(s) stands for the hard scattering 
improved Born cross-section at the reduced center-of-mass energy squared (see || and 
references therein). It should be emphasized that a(s) contributes like a distribution to 
c(s), since a(s) is smeared by the flux function over a wide range of the kinematic variable 
((f)(z;s) is peaked at z — 1 and has a tail till a cut at lower values of z). Expression (33) 
represents a(s) in a form with explicit integration over the virtualities of unstable particles 
(over the invariant masses of the corresponding final states). The quantity &o(s) reads 

- / \ fermion-loop-schcme / \ 

-a° n - shell > tIee ( s ; M+, Af_) x 2a ImS 2 (0)/ImS 1 (0) x TT S(s K - M 2 K ) x BR^ ee 

+ .on-shell,tree^ . M+? M ^ x j-j _ ^ x BR toee/boson-one-loop+real-photon 

K=± 

+ .on-shell, boson-one-loop + real-photon ^ . ^ ^ x n ^ _ ^ y ^tae _ (34) 

K=± 

Here the first term represents the result of the conventional fermion-loop scheme. All other 
terms describe its generalization. Factor 2 in the second term reflects the presence of two 
intermediate unstable particles of equal masses; BR± ee denotes their on-shell branching in the 
tree approximation. In the third term one BR is determined with the bosonic corrections to 
the partial width. (In fact the third term includes the sum of two subterms, with the modified 
BR for one of two unstable particle.) In (34) we have used the relation almSi(O) = MTq(M), 
with r (M) being the total on-shell width at tree level. This relation follows from unitarity 
and may be verified by direct calculation p2| . 



5 In general case, 4>(z; s) includes experimental cuts and hardware factors. 
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Formula (34), in principle, may be further simplified by carrying out the complete AO 
expansion in its first term. The key formula for this simplification, written in an OMS-like 
scheme of UV renormalization, is as follows: 



Wi(a; r) = [a ImS^O)] -1 7r<5(r) + Vp\ + 0(a). (35) 
Substituting (35) into (34) we finally get 

a (s ; s+, s-) = VP a on/off - she11 ' tree ( s ; M+, s-) x 5{s + -Ml) x BR? e 
+VP ff off/on - sM ' tree (s ; s+ , M_) x 5(s. -Ml ) x BRl™ 
-cr° n - shell ' tree ( S ;M + ,M_) x 2aImS 2 (0)/ImS 1 (0) x T7 S(s K - M 2 ) x BR? e 
+cr° n " shell ' tree (s;M + ,M_) x J] - M l) x ER K tree / boson - one - lo °P +real -P hoton 

K=± 

+5 .on-shell, tree + one-loop + real-photon^ . M+? M _) x J] ^ _ M ^ x ^tree _ (gg) 

Here the first two terms and the leading-order tree contribution in the last term originate 
from the fermion- loop-scheme term in formula (34). Namely, the first two terms accumulate 
the contributions originating from CC03 subprocesses and simultaneously from subprocesses 
of single vector boson production. In the case of CC03 subprocesses the symbol VP means 
that the corresponding off-shell unstable propagator squared is approximated, in accordance 
with (35), by VP 1/r 2 . Another unstable particle in this case is considered as real, produced 
on shell. In the case of single production the vector (unstable) boson is considered as real, 
too. The symbol VP in this case is superfluous and has to be omitted. 

It should be noted that in the above discussion the usage of prescription of the principal 
value in the case of the CC03 subprocesses is necessary, because the usage of any other 
prescription for 1/r 2 in (31) and (35) may change the results in (36). This remark, however, 
does not concern the anomalous term, since it arises from the contribution singular in a to 
Wu(a;r) (see (20) and the note after (22)). 

The above results may be easily generalized to other processes with unstable particles 
production, including the processes of the "NC"-type and of the "mixed" type. In the latter 
case there appear the additional terms in (34) and (36) of almost the same structure, that 
describe the ZZ-contributions (of course, configurations of Fig. ^ are not relevant in the 
ZZ-case). The generalization for the case of the differential cross-sections should not lead to 
difficulties either. However, the discussion of this question is beyond the scope of the present 
paper. 



8 Discussion 

In this paper we have found a practical method that ensures both the gauge cancellations 
and a fixed precision of description of processes mediated by the production and decay of 
unstable particles in electroweak theory. The solution has been found in the framework of the 
modified PT which consists of an expansion in the coupling constant directly of probabilities 
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instead of amplitudes In the general case and within any fixed precision the proposed 
method is based on the results obtained earlier in the framework of the background-field 
formalism. Within NLO precision and in the case of W- and/or Z-pair production we have 
found a way also in the usual formalism, basing ourselves on results of the fermion-loop 
scheme ||. But in contrast to the fermion-loop scheme we perform the description explicit 
taking into account the bosonic corrections and the two-loop corrections to the self-energy of 
unstable particles in the resonance region. From the practical point of view the latter result 
apparently is the most important one obtained in this paper. 

It should be noted that the result on gauge cancellations in the completely expanded 



probability generally was expected since the appearance of [13|. This work demonstrated 



that calculation of the probability of a process mediated by unstable particle production may 
be reduced to the regular fixed-order calculation. The latter fact gave a reason to think that 



the problem of the gauge invariance was solved, too (T3(]. However, on closer examination 
of the problem a large distance to obtaining the result was found. Indeed, the Feynman 
rules in the modified PT and in the conventional PT coincide only out of the mass shell. 
So, the gauge invariance a priori occurs only in this area of kinematic variable, but not 
in a neighborhood of the mass shell, where one has the non-standard contributions of the 
modified PT. However, due to the delta-functions these non-standard contributions do not 
vanish irrespectively of that as far as the neighborhood is small. Moreover, the properties 
of these non-standard contributions are unknown in advance. In particular, it is not known 
whether they are gauge invariant or not. This becomes especially clear in the higher orders 
of the AO expansion, beginning with NLO, where the non-standard contributions include 
the derivatives of the delta-functions. 

In the present paper we avoid this difficulty by making use of incomplete AO expansions 
of the unstable propagators squared. This method allows us to carry out a comparison of 
the results of the modified PT with those of the conventional PT with Dyson resummation. 
The special role in this scheme is assigned to the anomalous additive term that corrects the 
result of incomplete expansions of the propagator squared in the vicinity of mass shell. 

We draw attention to a particular problem that arises when comparing the results ob- 
tained in the modified PT with the results obtained in the conventional approach with use 
of the background-field formalism. (It should be emphasized that there is no such problem 
in the conventional formalism on using the fermion-loop scheme.) This is the problem of the 
residual dependence on the quantum gauge parameter which still remains in the background- 



field formalism after all gauge cancellations [ rOfl . In the general case this dependence may 
pass on to the results of the modified PT. Nevertheless, in the probability owing to the phe- 
nomenon of changing of the order of particular contributions, which in turn is a corollary of 
the singularity in the coupling constant in the propagator squared, one may expect that this 
dependence will drop out within the given precision. Note that we indicate only an oppor- 
tunity of giving a solution to the problem, but not the actual solution. However, earlier it 
was not clear even how to initiate a solution || 10|, 12 |. 

Now let us discuss the generalization of the fermion-loop scheme. Remember that it 
provides both the NLO precision and the gauge invariance. As a matter of fact the gener- 
alization consists in the instruction on how to use the already known results of the former 
calculations. Indeed, the first term in formula (34) and the last two terms are already known 
The remaining second term is actually known too, since ImS 2 (0) in it may be re- 
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placed by a one-loop correction to the width of W-boson divided by its mass (see Property 3 
and discussion thereof). 

The last two terms in (34) present exhaustive description of the bosonic corrections, which 
formerly were a "stumbling block" for gauge invariance. In the proposed generalization they 
are located in the ^-matrices squared, with the one S-matrix is for the on-shell production of 
unstable particles and the other ones are for their decays. The key observation that leads us 
to this simple result is the absence in the NLO approximation of non-factorizable corrections 
of Fig. H (see Sect. 0). This property demonstrates one of the differences of the approach 
offered with the well known double-pole approximation (DPA) [jj. 

In fact there are more serious differences, too. The most obvious one is as follows. By 
definition, DPA does not take into account the single-pole contributions, while the fermion- 
loop scheme, and consequently our generalization, does. The mentioned missed contributions 
are of order T\y/Mw, which amounts to 2-3%. This is too much compared with the current 
accuracy of LEP2 results. In view of this problem, DPA usually is applied for the calculation 
of the radiative correction only, but not of the Born term. The Born term is considered 
in the framework of the conventional PT with "by-hand" substitution of Breit-Wigner's 
propagators instead of the free propagators for unstable particles [0]. However, the latter 
operation violates the gauge invariance and leads to some error, too, and a certain estimate 
of this error does not exist. All that it is possible to say is that this error most likely is 
somewhat less than the shift of the amplitude caused by this substitution, but it is not 
clear how far less. Nevertheless, the latter quantity again is 0(Tw/Mw)- So, the discussed 
approximation of the Born term cannot be considered satisfactory. 

Another serious difficulty of DPA is its inapplicability in the vicinity of phase-space 
boundaries. This effect arises as a result of the "mapping" procedure in the phase space or 
the "analytic continuation" of the constant residues at the double pole in the amplitude 
First of all, it manifests itself near the threshold region, where the DPA uncertainties are 
blowing up. Another evident restriction of applicability of DPA is the region lying far off from 
resonance where the pole-scheme expansion cannot be viewed as an effective expansion in 
powers of T^/Mt^. Remember that the proposed generalization of the fermion-loop scheme 
does not need the mapping procedure and does not use the pole expansion. Consequently, 
it should be free from these difficulties. 

In conclusion, let us discuss two problems that arise in the modified PT approach. In 
fact, the topic is about a comparison of the results of the modified PT with Breit-Wigner's 
parameterization of unstable particles. The first problem concerns the definition of the 
"physical" mass and the width of unstable particles. In this connection it should be noted, 
first of all, that both these quantities are pseudo-observables that are to be determined on 
the basis of realistic observables 0, p2] . Next, let us note that in the framework of the 
modified PT both these quantities are secondary ones, which should be determined on the 
basis of primary objects, such as the renormalized Lagrangian mass, coupling constant, etc. 
Apparently, the most radical way consists of the identification of both these quantities with 
the position of the pole in the complex plane of the full unstable propagator (in the spirit 
of the pole scheme). This procedure is equivalent to finding a solution to the equation 
M 2 — s p — £(%) = with s p = Mp—iMpTp. This operation may be done perturbatively 



The second problem concerns the presence of the delta-functions and VP's in the results 
of the complete AO expansion. (The actuality of this problem is not so high while making use 
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of incomplete AO expansions as employed in this paper.) The problem may be designated 
as an illusory discrepancy between the presence of these singular functions in the amplitude 
squared and the notion of a continuity of the physical results as functions of the kinematic 



variables. This problem has been stated and discussed in [13| . The essence of the solution is 
reduced to the observation that actually there is not a direct identity between the squared 
amplitude (that is usually calculated) and the genuine probability (that is observed). In fact, 
a formal expression for the probability, calculated on the basis of the unitarity diagrams, 
should necessarily be subjected to the operation of integration, or "smearing" with some 
weight function before it becomes an observable quantity^ In processes with light and 
charged initial and/or final states such a smearing function is, first of all, the flux function 
that extracts and exponentiates the infrared and collinear singularities from the cross-section 
(see formula (32)). 

In the case of W- and/or Z-production the smearing is large enough: on some distance 
from threshold its final effect is about ten or more percent and it covers an area which is 
certainly greater than the unstable-particle width (see formulas (32) and (33)). Consequently, 
after the convolution the contributions of the delta-functions and VP's become strongly 
distributed over a wide area. Mathematically, this fact means that the smearing may be 
considered as an effect of order 0(1), which is required in the AO. Therefore, even in the 
case of ideal energy resolution in a given experiment, the presence of the singular functions 
in the results of the AO expansion becomes invisible, at least within the given precision 
of the description. The above conclusion should hold both for the total and differential 
cross-sections, since the differential cross-sections are determined through the convolution 
procedure, too |]J. 
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